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ABSTRACT 

We show that the star known as 2MASS J05162881+2607387 (hereafter J0516) is a double- 
lined eclipsing binary with nearly identical low-mass components. The spectroscopic elements 
derived from 18 spectra obtained with the High Resolution Spectrograph on the Hobby-Eberly 
Telescope during the Fall of 2005 are K 1 = 88.45 ± 0.48 km s" 1 and K 2 = 90.43 ± 0.60 km s"\ 
resulting in a mass ratio of q — K1/K2 — 0.978 ±0.018 and minimum masses of Mi sin i = 
0.775 ± 0.016 M Q and M 2 sin 3 i = 0.759 ± 0.012 M Q , respectively. We have extensive differential 
photometry of J0516 obtained over several nights between 2004 January-March (epoch 1) and 
2004 October-2005 January plus 2006 January (epoch 2) using the lm telescope at the Mount 
Laguna Observatory. The source was roughly 0.1 mag brighter in all three bandpasses during 
epoch 1 when compared to epoch 2. Also, phased light curves from epoch 1 show considerable out- 
of-eclipse variability, presumably due to bright spots on one or both stars. In contrast, the phased 
light curves from epoch 2 show little out-of-eclipse variability. The light curves from epoch 2 and 
the radial velocity curves were analyzed using our ELC code with updated model atmospheres 
for low-mass stars. We hnd the following: Mi = 0.787 ± 0.012 Mq, R 1 = 0.788 ± 0.015 R Q , 
Mi = 0.770 ± 0.009 M Q , and R 2 = 0.817 ± 0.010 R Q . The stars in J0516 have radii that are 
significantly larger than model predictions for their masses, similar to what is seen in a handful 
of other well-studied low-mass double-lined eclipsing binaries. We compiled all recent mass and 
radius determinations from low-mass binaries and determine an empirical mass-radius relation 
of the form R(R Q ) = 0.0324 + 0.9343M(M Q ) + 0.0374M 2 (M Q ). 

Subject headings: binary systems: low mass stars, individual (2MASS J05162881+2607387) 



1. Introduction 

Binary stars offer the best opportunity for ac- 
curate measurements of the radii and masses of 
stars, measurements essential to verify stellar evo- 
lution theory and to determine the properties 



1 Based on observations obtained with the Hobby-Eberly 
Telescope, which is a joint project of the University of 
Texas at Austin, the Pennsylvania State University, Stan- 
ford University, Ludwig-Maximilians-Universitat Miinchen, 
and Georg-August-Universitat Gttingen. 



of other diverse objects such as white dwarfs, 
neutron stars, black holes, and extra-solar plan- 
ets. These measurements rely on dynamical con- 
straints from observed radial and rotational veloc- 
ities, and geometric constraints from photometric 
time series observations. Using computer mod- 
els of binary stars we can derive physical parame- 
ters (e.g. masses, radii, etc.). Understanding the 
structure and evolution of stars is a basic goal of 
stellar astronomy, and is required in most other 
branches of astronomy. Critical tests of evolution 
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Table 1: Mass and Radius of sample double- lined eclipsing binaries plotted in Figure 10. 



Name: 



Mass (M Q ) 



Radius (Rq 



Reference 



V818 Tan B 
RX J0239. 1-1028 
RX J0239. 1-1028 
GU Boo A 
GU Boo B 
YY Gem AB 
NSVS0103 A 
NSVS0103 B 
TrES-HerO-07621 
TrES-HcrO-07621 
BW3 V38 A 
BW3 V38 B 
CU Cnc A 
CU Cnc B 
CM Dra A 
CM Dra B 



0.7605 ±0.062 
A 0.736 ±0.009 
B 0.695 ±0.006 
0.610 ±0.007 
0.599 ± 0.006 
0.5992 ± 0.0047 
0.540 ± 0.002 
0.498 ± 0.002 
A 0.493±0.003 
B 0.489±0.003 
0.44±0.07 
0.41±0.09 
0.4333 ±0.0017 
0.3890 ±0.0014 
0.2307±0.0010 
0.2136±0.0010 



0.768 ±0.010 
0.735 ±0.018 
0.710 ±0.016 
0.623 ±0.016 
0.620 ±0.020 
0.6191 ±0.0057 
0.527 ±0.002 
0.505 ±0.002 
0.453±0.060 
0.452±0.050 
0.51±0.04 
0.44±0.06 
0.4317 ±0.0052 
0.3908 ± 0.0094 
0.2516±0.0020 
0.2347±0.0019 



Torres & Ribas (2002) 
Ribas (2005) 
Ribas (2005) 

Lopez-Morales & Ribas (2005) 
Lopez-Morales & Ribas (2005) 
Torres & Ribas (2002) 
Lopez-Morales et al. (in prep) 
Lopez-Morales et al. (in prep) 
Creevey et al. (2005) 
Creevey et al. (2005) 
Maceroni & Montalban (2004) 
Maceroni & Montalban (2004) 
Ribas (2003) 
Ribas (2003) 

Lacy (1977) and Metcalfe et al. 
Lacy (1977) and Metcalfe et al. 



(1996) 
(1996) 



theory for stars other than the Sun can be made 
on a small set of eclipsing binary stars [see, Pols 
et al. 1997; Schroder et al. 1997]. In general, when 
accurate tests are available, the results of stellar 
evolution models compare favorably to data for 
main sequence stars with masses greater than one 
solar mass (Pols et al. 1997). In contrast, evi- 
dence has been growing that the models for stars 
on the lower main sequence have problems when 
confronted with precise data from eclipsing bina- 
ries. For example, Torres & Ribas (2002) showed 
that all available evolutionary models underesti- 
mate the radii of the components of the M-star 
YY Gem by about 20% and overestimate the ef- 
fective temperatures by 150 K or more. Similar 
discrepencies are found in V818 Tau (Torres & 
Ribas 2002), in CU Cnc (Ribas 2003), and GU Boo 
(Lopez-Morales & Ribas 2005). See Table 1 for a 
compilation of recent mass and radius determina- 
tions derived from low-mass eclipsing binaries. 

The disagreement between the models and the 
data for these binaries is troubling since models 
for low mass stars are used to estimate the ages 
for open clusters and individual T Tauri stars by 
placing them in an HR diagram. Since the number 
of well-studied low-mass binaries is still relatively 
small, observations of additional low mass bina- 
ries would be extremely useful. The star known 
as 2MASS J05162881±2607387 (hereafter J0516) 
was discovered to be an eclipsing binary by Schuh 



et al. (2003), who noticed its variability during 
the course of an extensive campaign to monitor a 
pulsating white dwarf. The eclipse period was de- 
termined to be 1.29395(25) days. On the basis of 
a moderate resolution spectrum and photometric 
colors, Schuh et al. (2003) determined an effec- 
tive temperature of 4200 K for the primary, corre- 
sponding to a spectral type of K7V. Assuming that 
the orbital period is the same as the eclipse period, 
Schuh et al. (2003) modelled the light curve and 
derived a mass ratio of M 2 /Mi « 0.11, a radius 
ratio of R2/R1 ~ 1, and a temperature ratio of 
T2/T1 « 0.6. Based on their modelling results, 
Schuh et al. (2003) suggested that J0516 consists 
of a late K-type star paired with a brown dwarf. 
Schuh et al. (2003) also briefly considered the pos- 
sibility that the orbital period of J0156 was twice 
the eclipse period (i.e. 2x 1.29395 = 2.58790 days). 
In that case, the primary and secondary eclipses 
are identical in their phased light curve. Our spec- 
troscopic results clearly indicate that the orbital 
period of J0516 is in fact 2.58791 days. 

We discuss below our observations of J0516, 
models of the light and velocity curves, and im- 
plications of our results. 
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Riglit Ascension (J2000) 



Fig. 1. — A 10 x 10 arcminute /-band image of the field of J0516 showing the seven comparison stars (Table 
2). The position of J0516 is indicated by the hash marks. 



Table 2: USNO-B1 Catalouge Coordinates of the 
seven comparison stars. The star number indi- 
cates the star in Figure 1 



Star 


USNO Number 


a 


S 




(U1125) 


(J2000) 


(J2000) 


1 


02245259 


05 16 38.302 


26 11 34.98 


2 


02244488 


05 16 35.311 


26 03 18.40 


3 


02244635 


05 16 35.933 


26 05 44.20 


4 


02243339 


05 16 30.797 


26 08 09.85 


5 


02243426 


05 16 31.133 


26 09 33.66 


6 


02245357 


05 16 38.688 


26 09 32.18 


7 


02248141 


05 16 49.409 


26 08 03.88 



2. Observations and Reductions 
2.1. Differential Photometry 

J0516 was observed extensively over 44 nights 
from January 2004 to January 2006 at the Mount 



Laguna Observatory (MLO, located in Southern 
California). We used the lm telescope equipped 
with a Loral 2048 x 2048 CCD (binned 2x2, yield- 
ing a pixel scale of 0.8 arcseconds per pixel) and 
standard V, R and I filters. Typically, the integra- 
tion times were 300-360 seconds, 240-300 seconds, 
and 180 seconds for V, R 1 and I bandpasses, re- 
spectively, with each observing session lasting sev- 
eral hours per night. These observations were di- 
vided into two main epochs based on the source 
behavior, January - March 2004 (epoch 1) and 
October 2004 - January 2005, plus January 2006 
(epoch 2). In epoch 1, there are 398 I band im- 
ages, 369 R band images, and 137 V band images. 
There were fewer V band images because of the 
difficulty of obtaining quality images through the 
occasional clouds and periods of poor seeing. In 
the epoch 2, there are 678 / band images, 588 R 
band images, and 481 V band images. 
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0.5 1 1.5 2 

Phase (P = 2. 58791 d) 

Fig. 2. — The light curves from the two epochs. The star was more active in epoch 1 when compared to the 
"quiesent" period seen in epoch 2, namely the star was approximately 0.1 magnitude brighter in all three 
bandpasses. 



Standard IRAF 2 tasks were used to remove the 
electronic bias and to perform the flat-fielding cor- 
rections. The IRAF task 'imalign' was used to re- 
move the differences in the pixel locations of the 
stellar images and to place all the CCD images on 
the same relative coordinate system. 

Differential light curves were derived using a 
two step process. Stable comparison stars were 
found using an automated script that utilizes Stet- 
son's programs DAOPHOT HE, ALLSTAR, and 
DAOMASTER (Stetson 1987; Stetson, Davis, & 
Crabtree 1991; Stetson 1992a & 1992b). To do 
this, the point spread function (PSF) was deter- 
mined from each image from fits to several iso- 

2 IRAF is distibuted by the National Optical Astronomy Ob- 
servatory, which are operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under the coopera- 
tive agreement with the National Science Foundation 



lated stars, and the instrumental magnitudes were 
computed using an aperture with a radius of 6 
pixels. We identified seven stable stars that are 
relatively bright and isolated. The coordinates of 
these seven stars are given in Table 2 and a finding 
chart is shown in Figure 1. Then, the IRAF task 
'phot' was used to perform aperture photometry 
on J0516 and the seven comparison stars using a 
sequence of 14 concentric apertures of radius 3- 
16 pixels. The IRAF implimentation of Stetson's 
curve-of-growth technique (Stetson 1990) was used 
to derive optimal instrumental magnitudes cor- 
responding to the largest aperture. Finally, the 
light curve of J0516 was found differentially by 
using the average instrumental magnitudes of the 
seven comparison stars. The phased differential 
light curves from both epochs are shown in Fig- 
ure 2. The mean magnitudes of J0516 deter- 
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mined by Schuh et al. (2003) are V = 18.1 ± 0.1, 
R = 16.8 ±0.3, and I = 15.84 ±0.3. We note that 
since the source is very red, the signal-to-noise is 
highest in I and lowest in V. 

2.2. Echelle Spectroscopy 

We have obtained 18 spectra (R — 15, 000) 
using the High Resolution Spectrograph (HRS, 
Tull 1998) and the Hobby-Eberly Telescope (HET, 
Ramsey et al. 1998). The instrument configu- 
ration consisted of a resolving power of 15,000, 
the central echelle rotation angle, the 316 groove 
mm -1 cross disperser set to give a central wave- 
length of 6948 A, the 2 arcsecond science fiber, 
and two sky fibers. The source was observed by 
the HET staff in dark skies and good seeing con- 
ditions between 2005 October 10 and December 
27. The exposure times were 1800 seconds each, 
split up into two 900 second parts to aid in the re- 
moval of cosmic rays. We also have the spectrum 
of the bright K7V star HD 28343 for comparison 
purposes, as well as ten spectra of six radial veloc- 
ity standard stars selected by the HET staff from 
the list of Nidever et al. (2002). 

The electronic bias was removed from each im- 
age. Then the pairs of 900 second exposures were 
combined using the 'crreject' option for cosmic ray 
removal. This worked very well, and the resulting 
spectra were largely free from cosmic rays. The 
echelle spectra were extracted and wavelength cal- 
ibrated using the IRAF 'echelle' package. The 
signal-to-noise ratios of the extracted spectra gen- 
erally lie in the range of about 10 to 20 per pixel. 
The HRS detector is a mosiac of two CCDs, and 
owing to the severe fringing on the "red" CCD, 
we have utilized only the portion of the spectra 
imaged on the the "blue" CCD, where the wave- 
length coverage is about 5100-6900 Angstroms. 

The stability of the instrument was checked by 
cross correlating the spectra of the standard stars 
against each other using the 'fxcor' task in IRAF. 
The standard deviation of the differences between 
the velocities measured using fxcor and the veloc- 
ities given in Nidever et al. (2002) was about 1.7 
km s^ 1 . 

Radial velocities for the J0516 spectra were de- 
rived using the "Broadening Function" technique 
developed by Rucinski (1999). This technique for 
extracting radial velocities of double-lined bina- 
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Fig. 3. — Representative broadening functions 
(BFs) for four observations are shown with filled 
circles. Double Gaussian fits to the broadening 
functions are shown with the solid lines. The cross 
correlation functions (CCFs) are shown with the 
dashed lines. The orbital phases are indicated in 
the upper right of each panel, and the vertical dot- 
ted lines denote the systemic velocity of the binary 
(17.1 km s _1 ). Note how the BFs have thinner and 
more separated peaks than the CCFs. 

ries is often more robust than a simple cross cor- 
relation technique (Tonry & Davis 1979), espe- 
cially when the velocity separation between the 
components is on the order of the spectral res- 
olution. The spectrum of HD 28343 was used 
as the high signal-to-noise, sharp-lined reference 
spectrum. In preparation for the analysis, each 
echelle order was normalized to its continuum level 
using a three piece cubic spline. The normalized 
spectra were Doppler corrected to the heliocen- 
tric rest frame and the echelle orders were merged 
using a linear dispersion in the wavelength range 
5898.68 < A < 6689.00 A with a pixel size of 
0.1168 A. In fifteen out of the eighteen spectra, the 
broadening functions (BFs) have two clear peaks, 
indicating that the spectrum is double-lined. The 
remaining three spectra happened to be taken at 
conjunction phases, and hence resembled single- 
lined spectra. Radial velocities were measured us- 
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Fig. 4. — The radial velocity curve of the two 
stars. 

ing Gaussian fits to the peaks in the BFs. Figure 3 
shows representative BFs and the double Gaussian 
fits, and the corresponding cross correlation func- 
tions (CCFs) for comparison, where the spectrum 
of HD 28343 was used as the template spectrum. 
The peaks in the BFs are much better separated 
than the peaks in the CCFs, resulting in more re- 
liable radial velocities. Figure 4 shows the radial 
velocities for both components folded on the pe- 
riod determined below. 

The full widths at half maximum of the Gaus- 
sian fits to the J0516 broadening function peaks 
average about 16 km s _1 , compared to an aver- 
age of 11 km s _1 for the resulting peaks for the 
various sharp- lined radial velocity standard stars. 
This indicates that the rotational velocities of the 
stars in J0516 are marginally resolved (e.g. the re- 
solving power corresponds to a velocity of w 20 
km s _1 ), and are on the order of 16 km s^ 1 . 

We constructed "restframe" spectra for both 
stars by Doppler correcting each of the 18 individ- 
ual spectra to zero velocity and averaging them 
together using 'minmax' rejection. This was done 
for the primary, which we define as the compo- 
nent that is eclipsed at HJD 2452251.5173 [the T 
given in Schuh ct al. (2003)] and then again for 
the secondary. The lines of the secondary appear 
to be removed from the restframe spectrum of the 
primary (and vice versa) reasonably well using the 
minmax rejection. Note, however, that the lines 
in a given restframe spectrum are diluted by the 
continuum from the other component. Figure 5 
shows six different cchellc orders of each restframe 



spectrum. Each order has been normalized to its 
continuum, smoothed using a running average of 
five pixels, and the spectrum of the secondary has 
been offset by 0.25 units for clarity. For the sake 
of the presentation, the J0516 spectra showing the 
Mg I b feature near 5180 A have not been sky sub- 
tracted owing to the poor signal-to-noise in the 
object there. Also shown is the spectrum of the 
K7V star HD 28343 (obtained with the same in- 
strumental configuration), which has been normal- 
ized to its continuum and scaled by 0.5 in order 
to account for the dilution of the lines in the rest- 
frame spectra. Although the signal-to-noise ratios 
in the combined "restframe" spectra are not ter- 
ribly large, the two stars appear to have nearly 
identical line features, many of which resemble the 
lines seen in the K7V comparison. Thus the spec- 
tral type of K7V determined by Schuh ct al. (2003) 
appears to be correct. The relative areas under the 
BFs give the "luminosity" ratio, L 2 /Li, which is 
unity to within the errors. 

There is, however, very noteable difference be- 
tween the spectra of the J0516 components and 
the template spectrum: the stars in J0516 have 
no significant feature at Ha, either in emission 
or in absorption. In this regard, J0516 resembles 
some RS Canum Venaticorum and BY Draconis 
binaries that have "filled in" Ha line profiles (e.g. 
Fernandez-Figueroa, et al. 1994). Our photomet- 
ric light curves, discussed below, also indicate that 
one or both stars in J0516 have relatively high lev- 
els of stellar activity. The equivalent width of the 
Ha absorption line in the K7V template star is 
about 0.55 A, so presumably each star in J0516 
would have to have an Ha emission line with an 
equivalent width of about 0.55 A in order to pro- 
duce a roughly featureless spectrum at Ha. 

We examined the spectra near 6708 A for an 
indication of the Li I doublet, which is some- 
times used as an age indicator (e.g. Boesgaard & 
Tripicco 1986; Boesgaard & Budge 1988). There 
are no significant lines near this wavelength in ei- 
ther the primary restframe spectrum or the sec- 
ondary restframe spectrum, down to an equivalent 
width of about 0.02 A, which is roughly the noise 
level (see Figure 5). 

There are more sophisticated deblending tech- 
niques to analyze the spectra of double-lined bi- 
naries, but these typically require a grid of tem- 
plate stars. Thus, we are limited since the only 
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Fig. 5. — The six panels show the smoothed, continuum normalized spectra of the J0516 components and 
of the K7V star HD 28343. In each case, the spectra from top to bottom are the J0156 primary, the J0516 
secondary, offset by 0.25 units, and HD 28343, scaled by 0.5 to account for the dilution of the lines in the 
J0516 spectra. The J0516 spectra showing the Mg I b feature near 5180 A have not been sky-subtracted 
owing to the poor signal-to-noise there. 



appropriate template spectrum we have is that of 
the K7V star HD 28343 (the radial velocity stan- 
dards observed for us are all various subtypes of 
spectral class G). Nevertheless, we did construct 
an alternate restframe spectrum of the primary 
by scaling the normalized template spectrum by 
0.5, Doppler shifting it to match the velocity of 
the secondary for each observation, subtracting it 
from each observation, Doppler shifting the differ- 
ence spectra to remove the velocity of the primary, 
and finally, averaging the "difference spectra" to 
produce the final result. Apart from having ab- 
sorption lines a factor of two stronger, this alter- 
nate restframe spectrum looked very similar to the 



spectrum shown in Figure 5. The noise levels in 
the continua are similar, and there is no feature at 
Ha. 

3. Analysis 

3.1. System Variablity 

The light curves from epoch 1 are distinctly 
different than the light curves from epoch 2. As 
shown in Figure 2, the source was roughly 0.1 mag 
brighter in all three bandpasses during epoch 1 
when compared to epoch 2. Although we do not 
have standard star observations, an inspection of 
the instrumental V — I colors indicates that the 
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Phase (P = 2. 58791 d) 

Fig. 6. — The epoch 2 data with the best-fitting model from ELC and the corresponding residuals. Left: the 
epoch 2 spotless model and residuals. Right: the epoch 2 spot model with a single spot on each star. The 
bandpasses are, from top to bottom, /, R, and V. 



system was « 0.02 mag redder in V — I during 
epoch 1. Also, phased light curves from epoch 

1 show considerable out-of-eclipse variability, pre- 
sumably due to bright spots on one or both stars. 
In contrast, the phased light curves from epoch 

2 show little out-of-eclipse variability. Asymme- 
tries in the light curves seem to be a common fea- 
ture of low-mass binaries (e.g. the GU Boo light 
curves shown by Lopez-Morales & Ribas 2005). 
However, the activity in J0516 seems to be at a 
much higher level compared to other well-studied 
low mass binaries, especially the w 10% increase 
on the system brightness observed in epoch 1. If 



this increase in brightness is due to bright spots, 
then one would need to have a substantial portion 
of the surface of one or both components covered 
with bright areas in order to have them visible at 
all phases. 

3.2. Light Curve Modeling 

We modeled the light curves using our ELC 
code (Orosz & Hauschildt 2000) with updated 
model atmospheres for low mass main sequence 
stars and brown dwarfs (Hauschildt, private com- 
munication). Our "base" model has nine free pa- 
rameters: the inclination i, the mass of the pri- 
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mary Mi , the if- velocity of the primary K\ , the 
radius of the primary R\ , the temperature of the 
primary T\, the ratio of the radii R1/R2, the ratio 
of the temperatures T2/T1, the orbital period P, 
and the time of primary eclipse T . Note that for a 
given orbital period P and inclination i, specifying 
the primary mass Mi and if-velocity K\ uniquely 
determines the ratio of masses q — M 2 /Mi and 
the orbital separation a. We have found that in 
cases such this one, it is more efficient to explore 
the Mi,Ki parameter space rather than the q,a 
parameter space since the observed radial velocity 
curves give K\ and a rough value of M\ directly. 
We adopt a primary temperature of 4200 ± 200 K 
(Schuh ct al. 2003). The gravity darkening expo- 
nents were set according to the mean stellar tem- 
peratures according to the results of Claret (2000). 
"Simple" reflection was used (see Wilson 1990), 
assuming bolometric albedos of 0.5 for each star. 
We assume a circular orbit with synchronous ro- 
tation for both stars. We note, finally, that limb 
darkening coefficients are not needed since we have 
model atmosphere specific intensities tabulated at 
99 emergent angles. 

The light curves of close binary stars are some- 
times asymmetric about the conjunction phases, 
and these asymmetries are often attributed to 
spots (either bright or dark) on one or both com- 
ponents. The spots in ELC are parameterized in 
the same way as in the Wilson-Devinney (1971) 
code. They are circular regions specified by four 
parameters: the "temperature factor" Tf, the 
"latitude" of the spot center, the "longitude" of 
the spot center, and the angular radius of the spot. 
Bright spots have Tf > 1 and dark spots have 
Tf < I- 

We began by modelling the Epoch 2 light curves 
since these light curves are better sampled and 
show very little out-of-eclipse variability. The 
best-fitting model (with no spots) was arrived at 
by iteration and brute force. Our "observables" 
are light curves in V, R, and /, radial velocity 
curves for both stars, and a "luminosity ratio" of 
L 2 /Li = 1.00 ± 0.03, taken to be in the V-b&nd. 
We ran ELC's genetic optimizer several times to 
arrive at the best intermediate solution using lib- 
eral ranges for the free parameters. Using this so- 
lution, we scaled the error bars on the photometry 
and radial velocities to give x 2 = — 1 for each 
data set separately. The median error bars after 




0.5 1 1.5 



Fig. 7. — /-band residuals from epoch 2 model fits 
are shown phased on the orbital period and binned 
into 50 bins. The residuals marked (a) are from 
the best-fitting model which had a single spot on 
each star. The residuals marked (b) are for the 
model with no spot, and the residuals marked (c) 
are for the model with a single spot on the primary 
(see Tables 3 and 4 for parameters). Note the 
persistent feature near the primary eclipse (phase 
0). 

this scaling were 0.030 mag for V, 0.018 mag for R, 
0.017 mag for /, 1.16 km s _1 for the primary radial 
velocity curve, and 1.15 km s _1 for the secondary 
radial velocity curve. The data sets with the scaled 
error bars were optimized again using the genetic 
optimizer and a simple "grid search" technique. 
The best-fitting spotless model is shown in the left 
panels of Figure 6. 

We used the procedure outlined in Orosz et al. 
(2002) to find approximate la, 2a, and 3a confi- 
dence intervals. To estimate uncertainties on fit- 
ted and derived parameters we projected the 9- 
dimcnsional x 2 function into each parameter of 
interest. The la, 2a, and 3c confidence limits 
taken to be the ranges of the parameter where 
X 2 < Xmin + l' +4> an d +9, respectively. Since the 
genetic ELC code samples parameter space near 
Xmin extensively, computing these limits is simple. 
ELC saves from every computed model the \ 2 of 
the fit, the value of the free parameters, and the 
astrophysical parameters (e.g., the primary star 
mass, the radii of the components, etc.). One can 
then select out the lowest x 2 at each value of the 
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Table 3: The results of ELC modeling. 



Parameter 


Epoch 2 without 


Epoch 2 with 


Epoch 1 with 




spot 


spot on each star 


with spot 


Period (days) 


2.58791 ±0.00001 


2.58792 ± 0.00001 


2.58791 ± 0.00002 


T (HJD+2,450,000) 


2251.512 ±0.005 


2251.508 ±0.005 


2251.517 ± 0.005 


Inclination (deg) 


84.1 ± 0.1 


84.3 ± 0.1 


84.08 ± 0.04 


T2/T1 


0.989 ± 0.005 


0.988 ± 0.005 


0.981 ± 0.004 


a (R Q ) 


9.19 ± 0.05 


9.19 ± 0.05 


9.24 ± 0.10 


log(ffi) 


4.51 ± 0.02 


4.54 ± 0.02 


4.58 ±0.03 


1 / \ 

l0g(ff2) 


4.49 ± 0.01 


4.50 ± 0.01 


4.51 ± 0.02 


V\ sin-i (km s _1 ) 


15.93 ±0.28 


15.34 ±0.30 


14.88 ±0.50 


V 2 sini (km s _1 ) 


15.95 ±0.19 


15.89 ±0.20 


15.88 ±0.25 


K x (km s- 1 ) 


88.42 ±0.48 


88.45 ±0.48 


88.30 ±0.40 


K 2 (km s- 1 ) 


90.46 ±0.64 


90.43 ±0.60 


90.30 ±0.60 


7 (km s _1 ) 


17.1 ± 1.7 


17.0 ± 1.7 


17.0 ±0.1 


q = M 2 /M 1 


0.978 ±0.019 


0.978 ±0.019 


0.978 ±0.019 


My (M Q ) 


0.787 ±0.012 


0.787 ±0.012 


0.786 ±0.010 


M 2 (M ) 


0.770 ± 0.009 


0.770 ± 0.009 


0.769 ±0.008 


i?i (/?©) 


0.818 ±0.015 


0.788 ±0.015 


0.769 ±0.015 


R 2 (-R©) 


0.820 ±0.010 


0.817 ±0.010 


0.806 ±0.010 



parameter of interest. The values and their un- 
certainties of various fitting and derived parame- 
ters for the best-fitting Epoch 2 spotless model are 
given in Table 3. 

We then checked for systematic trends by 
searching for periodicities in the /-band residu- 
als (the /-band light curve has the best signal- 
to- noise). The Lomb-Scargle periodogram (Lomb 
1976; Scargle 1982) showed significant power at 
the orbital period. Indeed, a small modulation 
with an amplitude of about 1% is evident in the 
/-band residuals when they are phased on the or- 
bital period and binned into 50 phase bins (see 
Figure 7). There is also a small systematic prob- 
lem with the model near primary eclipse. 

We performed additional fits to the Epoch 2 
light curves using models with spots. We had five 
different spot configurations: a spot on the pri- 
mary only, a spot on the secondary only, two spots 
on the primary only, two spots on the secondary 
only, and a single spot on both the primary and 
secondary. All of the models with spots had x 2 
values that are significantly lower than the spot- 
less model (see Table 4). The phased and binned 
/-band residuals for the spotted models look very 
similar to each other. Figure 7 shows the phased 
and binned residuals for two of the spot models. 



The low-level modulation near phase 0.5 is gone, 
but the feature near primary eclipse at phase re- 
mains. Apart from some small differences in the 
derived radii of the stars, the model and derived 
parameters for the spotted models agree quite well 
with the spotless model. As an example, Table 3 
gives the parameters for the model with a single 
spot on each star, and the panels on the right side 
of Figure 6 show this model. 

The temperature of the primary is not con- 
strained by our light curves alone. However, the 
ratio of the temperatures is quite well constrained 
by the light curves, and is quite close to unity (see 
Table 3), as we expected based on the spectro- 
scopic results (i.e. the restframe spectrum of the 
primary looks the same as the restframe spectrum 
of the secondary to within the noise level) . 

In all of the modelling described previously, we 
used the observed "luminosity ratio" in the V- 
band of 1.0 ± 0.03 as an additional constraint. 
Basically, we compute the quantity Xi urn = (1 — 
L 2 /L 1 ) 2 /0m 2 and add it to the total \ 2 , which 
has the effect of selecting models that have L 2 ~ 
L\. This is an entirely reasonable constraint since 
the spectroscopic results indicate that the primary 
star has the same spectral type as the secondary 
star, and that the two stars have nearly equal 
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0.5 1 1.5 2 

Phase (P = 2. 58791 d) 

Fig. 8. — The epoch 1 data with the best-fitting model from ELC and the corresponding residuals. See 
Table 3 for the parameters. The bandpasses are, from top to bottom, /, R, and V. 



masses (i.e. K\ « K2). In order to test for possi- 
ble systematic errors, we did additional model fits 
to the epoch 2 light curves without the luminos- 
ity ratio constraint. For the baseline model with 
no spots, we found basically the same tempera- 
ture ratio (T 2 /Ti = 0.984) as before. However, 
the ratio of the radii was R1/R2 = 1.064, which 
yields a l/-band luminosity ratio of L1/L2 ~ 1.5. 
It seems unlikely that two stars with nearly equal 
masses and temperatures should have such dissim- 
ilar radii, so we reject this model and only con- 
sider the models that have the luminosity ratio 
constraint imposed. 

The most interesting astrophysical parameters 
are the masses and radii of the stars, and Table 4 
gives these values for all of the Epoch 2 models. 
Formally, the model with a single spot on each 
star gives the best fit. While the derived masses 
depend mainly on the radial velocity curves, the 
light curves have some weak dependence on the 
scale of the binary since ELC uses model atmo- 
sphere intensities that are tabulated in tempera- 
tures and gravities. As a result, the derived masses 



for each of the models are not exactly the same. 
There is a maximum spread of 0.007 M Q for the 
mass of the primary and 0.006 M Q for the sec- 
ondary, respectively, which is is less than the for- 
mal la errors. Based on the best-fitting model, 
we adopt values of Mi = 0.787 ± 0.012 Mq and 
M 2 = 0.770 ± 0.009 M for the primary and sec- 
ondary masses, respectively. On the other hand, 
the derived radii show a bit more spread be- 
tween the various models. There is a maximum 
spread of 0.020 Rq for the radius of the primary 
and 0.013 Rq for the radius of the secondary, re- 
spectively. Based on the best-fitting model, we 
adopt values of R x = 0.788 ± 0.015 R & and R 2 = 
0.817 ± 0.010 Rq for the primary and secondary 
radii, respectively. We caution that these values 
for the radii probably have systematic errors on 
the order of 0.02 Rq, given the spread in radii for 
the various models, the persistent feature in the 
residuals near phase shown in Figure 7, and the 
need to impose the spectroscopically determined 
luminosity ratio. 

For completeness, we also fit the epoch 1 data 



11 



Table 4: Masses and radii for Epoch 2 models 



model 


Mi 


M 2 


Ri 


i?2 


x 2 


comment 


no spots 


0.787 ±0.012 


0.770 ±0.009 


0.818 ±0.015 


0.820 ±0.010 


1774 




1 spot each 


0.787 ±0.012 


0.770 ±0.009 


0.788 ±0.015 


0.817 ±0.010 


1551 


dark and bright spot 


2 spots, primary 


0.792 ±0.012 


0.774 ±0.009 


0.797 ±0.015 


0.819 ±0.010 


1570 


dark spots 


2 spots, secondary 


0.787 ±0.012 


0.769 ±0.009 


0.795 ±0.015 


0.823 ±0.010 


1576 


bright spots 


1 spot, primary 


0.787 ±0.012 


0.770 ±0.009 


0.797 ±0.015 


0.810 ±0.010 


1603 


dark spot 


1 spot, secondary 


0.785 ±0.012 


0.768 ±0.009 


0.814 ±0.015 


0.818 ±0.010 


1604 


dark spot 


epoch 1 


0.786 ±0.010 


0.769 ±0.015 


0.769 ±0.015 


0.806 ±0.010 




bright spot 



2000 3500 3000 3500 

Time (HJD-2450000) 

Fig. 9.— O-C diagram of J0516. The Xs indicate 
data from this work. The stars indicate data from 
Schuh ct al. (2003). 

using the same basic model as above and adding 
a single spot on the primary. The parameters 
are given in Table 3 and the phased light curves, 
the best-fitting models, and residuals are shown 
in Figure 8. The model fits reasonably well, and 
for the most part the fitted and derived param- 
eters agree quite well with those found from the 
epoch 2 fits, in spite of the fact that the secondary 
eclipse is not well sampled in epoch 1, and, as 
noted above, there is considerable out-of-eclipsc 
variability in epoch 1. The one noteable exception 
is the primary radius, which is 0.769 ± 0.015 i? Q 
from the epoch 1 model and 0.788 ±0.015 Rq from 
the epoch 2 model, a difference of nearly 0.02 i? Q . 
Although the epoch 1 residuals have more scat- 
ter than their epoch 2 counterparts, the epoch 1 
residuals show no significant power at the orbital 
period in a Lomb-Scargle periodogram. 



3.3. Eclipse Timings 

We determined eclipse times by fitting a 
parabola to the /-band observations near times 
of eclipse (the /-band observations have the high- 
est signal-to- noise) . The timings and cycle counts 
are given in Table 5 and have typical uncertain- 
ties of about 10 minutes. The table also gives the 
eclipse times determined by Schuh et al. (2003) 
with updated cycle counts. Schuh ct al. (2003) 
give no uncertainties in their eclipse timings, but 
from an inspection of their plotted light curves, it 
seems likely that their timings are more accurate 
than ours, so we assign their timing three times 
more weight than our timings. These timings are 
described by the linear ephemeris 

min I = HJD 2, 452, 251.518(1) + 2.587890(7)£ 

(error in the last digit is in parenthesis). This pe- 
riod differs by « 2a from the spectroscopic/photometric 
period derived from the ELC modelling. An O — C 
diagram is shown in Figure 9. The maximum de- 
viation is on the order of 20 minutes, and there is 
no obvious trend. 

4. Summary and Discussion 

We have derived fairly accurate values of the 
masses and radii of the two stars in J0516. The 
masses are determined to 1.5% precision and the 
radii are formally determined to 1.9% precision, 
although there may be a systematic error of a 
few percent on the radius determinations. This 
high level of precision compares favorably to the 
precision achieved on the small number of other 
low-mass eclipsing binaries mentioned in Section 
1. For low-mass stars with well determined masses 
and radii, the observed radii are consistently larger 
than what evolutionary models predict for their 
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Table 5: The HJD of the eclipse times, epochs, and 0-C values. 



Eclipse Time (HJD) Epoch 0-C (days) 



Reference 



2452251 


5164 


0.0 


-0.0013 


2452255 


4007 


1.5 


0.0012 


2452263 


1615 


4.5 


-0.0017 


2452264 


4599 


5.0 


0.0028 


2452265 


7487 


5.5 


-0.0024 


2452267 


0450 


6.0 


0.0000 


2452586 


6506 


129.5 


0.0012 


2452590 


5304 


131.0 


-0.0008 


2453018 


8301 


296.5 


0.0031 


2453022 


7100 


298.0 


0.0012 


2453048 


6001 


308.0 


0.0124 


2453287 


9700 


400.5 


0.0025 


2453304 


8000 


407.0 


0.0112 


2453357 


8401 


427.5 


-0.0004 


2453388 


8799 


439.5 


-0.0153 


2453753 


7800 


580.5 


-0.0076 



Schuh et al. (2003) 



This 



work 



masses. The stars in J0516 are no exception. Fig- 
ure 10 shows the stars in J0516 and those of other 
well-studied low-mass binaries in a mass-radius 
plot. We also show in Figure 10 theoretical mass- 
radius relationships with ages of 50 Myr, 100 Myr, 
and 300 Myr taken from the models of Baraffe 
ct al. (2002). No single model passes through 
all of the observations, and all of the stars are 
above the 300 Myr model. We formed an empiri- 
cal mass-radius relationship by fitting a parabola 
using weights on both the mass and radius through 
the other stars and the Sun (excluding J0516) and 
found the following equation 

R(R Q ) = O.O324+O.9343M(M Q )+O.O374M 2 (M ). 

Compared to this empirical relationship, the J0516 
primary agrees quite well, whereas the J0516 sec- 
ondary is slightly too large. 

Why is the secondary of J0516 so large for is 
mass? Ribas (2005) raises the possibility that the 
high level of stellar activity observed in the well- 
studied double-lined low-mass eclipsing binaries 
is what causes these stars to have radii that are 
larger than what is predicted. As we discussed 
earlier, J0516 is much more active than the other 
binaries, and the secondary's deviation from the 
mass-radius relationship is the most extreme. It 
is not clear, however, if higher levels of activity 
would automatically lead to larger radii. 

The Li I line is often used as an age indicator. 



An upper limit to the equivalent width of the Li I 
line at 6708 A of 20 mA translates roughly to an 
abundance of ^4(Li) = 0.0, which translates into a 
lower limit on the age of the binary of about 150 
Myr (Steinhauer 2003; Deliyannis private commu- 
nication). The evolutionary models of Baraffe et 
al. (2002) indicate that a 0.8 M star takes about 
70 Myr to contract and reach its "normal" main 
sequence radius. Taken at face value, the lower 
limit on the age of 150 Myr from the lack of lithium 
indicates that J0516 is well past its pre-main se- 
quence phase, and that the unusually large radius 
of the secondary cannot be explained by youth. 

The only practical way to establish the age 
of J0516 would be to associate it with a cluster 
or moving group. To do this, one would need 
good proper motions and distance for J0516 and 
a large number of stars in the nearby field. We 
know of no measurement of the proper motion of 
J0516 [for example, no proper motion measure- 
ment is given in the USNO-B catalog (Monet et 
al. 2003)]. The distance can be estimated from 
our model, but the result depends on the extinc- 
tion and on the assumed temperatures of the stars. 
Schuh et al. (2003) determined a color excess of 
E(B - V) = 0.9 ± 0.2 mag from their model 
of the spectral energy distribution. Assuming 
A v = 3.1E(B - V) and A K = 0.1UA V (Cardelli, 
Clayton, & Mathis 1989), the if -band extinction 
is 0.318 ± 0.017 mag. The absolute K magnitude 
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0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Mass (M ) 

Fig. 10. — Top: Mass versus radius of some observed "double-lined" eclipsing binaries mentioned in the 
text. J0516 is repesented by the two stars. Error bars smaller than the symbol size have been omitted. 
Evolutionary models from Baraffe et al. (2002) with ages of 50 Myr, 100 Myr, and 300 Myr are represented 
by the dot-dashed line, dotted line, and dashed lines, respectively. The solid line is weighted a quadratic fit 
to the observed data (including the Sun and excluding J0516). The secondary of J0516 has a larger radius 
than what the observed relationship defined by the other stars predicts. Bottom: The residuals from the 
quadratic fit. 



of the binary was found using filter-integrated sur- 
face brightnesses computed from NextGen mod- 
els by France Allard (private communication), as- 
suming the radii given in Table 3 and tempera- 
tures of 4200 ± 200 K. The average apparent mag- 
nitude in the if-band determined by the 2MASS 
survey is K = 13.113 ± 0.039 (see Schuh ct al. 
2003). We find a distance of d = 753 ± 34 pc. 
This distance is substantially larger than the dis- 
tances to the other well-studied low-mass binaries 
mentioned earlier [e.g. d = 140 ± 8 pc in the case 
of GU Boo (Lopez-Morales & Ribas 2005)]. 

J0516 is clearly a very important system in that 
it allows for further observational constraints for 
the lower main sequence. Owing to the enhanced 
variability, it would be worthwhile to obtain addi- 



tional photometry to better define the "quiescent" 
phase and to better establish the eclipse profiles in 
quiescence. In addition, moderate resolution but 
high signal-to-noise spectroscopy would be useful 
in order to further investigate the apparent lack 
of an Ha feature in either star and to check for 
possible correlations between the profile of the Ha 
feature (if any) and the activity level. 

The Hobby-Eberly Telescope (HET) is a joint 
project of the University of Texas at Austin, the 
Pennsylvania State University, Stanford Univer- 
sity, Ludwig-Maximilians-Universitat Munchcn, 
and Georg-August-Universitat Gottingen. The 
HET is named in honor of its principal bene- 
factors, William P. Hobby and Robert E. Eberly. 
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This research has also made use of the NASA/IPAC 
Infrared Science Archive, which is operated by the 
Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National 
Aeronautics and Space Administration, and also 
NASA's Astrophysics Data System. We thank 
Danielle Martino for the data taken at MLO from 
January - March 2004 and Con Deliyannis for 
helpful discussions. 
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